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The recent HESS detections of supernova remnant shells in TeV 7-rays con- 
firm the theoretical predictions that supernova remnants can operate as powerful 
<^ ! cosmic ray accelerators. If these objects are responsible for the bulk of galactic 

cosmic rays, then they should accelerate protons and nuclei to 10 15 eV and be- 
yond, i.e. act as cosmic PeVatrons. The model of diffusive shock acceleration 
allows, under certain conditions, acceleration of particles to such high energies 
and their gradual injection into the interstellar medium, mainly during the Sedov 
phase of the remnant evolution. The most energetic particles are released first, 
£j [ while particles of lower energies are more effectively confined in the shell, and 

are released at later epochs. Thus the spectrum of nonthermal paticles inside 
the shell extends to PeV energies only during a relatively short period of the 
evolution of the remnant. For this reason one may expect spectra of secondary 
7-rays and neutrinos extending to energies beyond 10 TeV only from T < 1000 yr 
old supernova remnants. On the other hand, if by a chance a massive gas cloud 
appears in the < 100 pc vicinity of the supernova remnant, "delayed" multi- TeV 
signals of 7-rays and neutrinos arise when the most energetic partices emerged 
t^. I from the supernova shell reach the cloud. The detection of such delayed emission 
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of multi- TeV 7-rays and neutrinos allows indirect identification of the supernova 
remnant particle PeVatron. 
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Subject headings: cosmic rays — gamma rays: theory — supernova remnants 
molecular clouds 



1. Introduction 

In 1934, Baade and Zwicky proposed that supernovae are responsible for the observed 
flux of cosmic rays (CR), provided that a substantial fraction of the kinetic energy released at 
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explosions of supernovae is converted into relativistic particles. Later, it has been recognised 
that relativistic particles can be effectively accelerated via Fermi mechanism at shock waves 
that form during the expan sion of supernova remnants (SNR) in the interstellar medium 



e.g. 



Malkov fc Drurvll200lh 



The particle acceleration in SNRs is accompanied by production of 7-rays and neutrinos 
due to interactions of accelerated CR protons and nuclei with ambient medi um. Recently, the 



HESS collaboration reported detection of some young SNRs in Te V T-ravs (lAharonian et al. 



2004h with fluxes quite close to the early theoretical predictions (IDrury. Aharonian fc Volki 



19941 ) . Though the HESS results constitute a very important advancements in the field, they 
still do not provide us with a definite and direct evidence of proton acceleration in SNRs; 
the competing inverse Compton scattering of directly accelerated electrons may significantly 
contribute to the observed 7-ray fluxes, provided that magnetic field in the acceleration 
region does not exceed 10 fiG. Production of tt° decay 7-rays is accompanied by radiation 
of neutrinos. Their detection, though challenging even for t he next generation of telescopes 
( ]Halzenll2007l ; iKappes et alll2007l ; ICostantini fc Vissanill2005l ) , would provide an unambiguous 
evidence for proton acceleration in SNRs. 

The observed CR spectrum extends without any distinct feature up to the so called 
knee at an energy of a few PeV. This suggests that the sources of galactic CRs, whichever 
they are, must be able to accelerate particles up to at least a few PeV. The standard lin- 
ear treatment of diff usive shock acceleration d oes not allow acceleration of protons beyond 
« 10 14 (-E>/3/iG) eV (ILagage Cesarskylll983l ) which falls short of the position of the knee, 
unless one assumes strong amplification of the magnetic field in the upstream region. Re- 
markably, the background magnetic field close to a shock which is efficiently ac celerating par- 
ticles can be amplified by a factor up to ~ 10 3 due to CR-driven instabilities iBell fc Lucek 
( 120011 ). The magnetic field amplification hypothe sis is supported also by observations of thin 



X-ray synchrotron filaments surroun ding SNR (IBamba et al.l 120031 ; IVink fc Laming! 12003 
Volk. Berezhko Ksenofontovll2005l ). 



A decisive and unambiguous indication of acceleration of PeV protons in SNRs can 
be provided by observations of 7-rays at energies up to 100 TeV and beyond. Because of 
the Klein-Nishina effect the efficiency of inverse Compton scattering in this energy band is 
dramatically reduced. Therefore unlike other energy intervals, the interpretation of 7-ray 
observations at these energies is free of confusion and reduces to the only possible mecha- 
nism - decay of secondary 7r°-meson. Although the potential of the current ground-based 
instruments for detection of such energetic 7-rays is limited, it is expected that the next 
generation arrays of imaging Cherenkov telescopes optimized in the multi-TeV energy range 
will become powerful tools for this kind of studies (lAharonian! 120051 ). 
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It should be noted that the number of SNRs currently bright in > 10 TeV 7-rays is ex- 
pected to be rather limited. Multi-PeV protons can be accelerated only during a relatively 
short period of the SNR evolution, namely, at the end of the free-expansion phase/beginning 
of the Sedov phase, when the shock velocity is high enough to allow sufficiently high ac- 
celeration rate. When the SNR enters the Sedov phase, the shock gradually slows down 
and correspondingly the maximum energy of the particles that can be confined within the 
SNR decreases. This determ ines the escape of the most energetic particles from the SNR 
Jptuskin & Zirakashvilj|20~0~5n . Thus, unless our theoretical understanding of particle accel- 
eration at SNR is completely wrong, we should expect an energy spectrum of CR inside the 
SNR approaching PeV energies only at the beginning of the Sedov phase, typically for a 
time < 1000 years. When the remnant enters the Sedov phase, the high energy cutoffs in 
the spectra of both protons and 7- rays gradually moves to lowe r energies, while the highest 
energy particles leave the remnant (IPtuskin fc Zirakashvilill2005l ). This can naturally explain 
why the 7-ray spectrum of the best studied SNR RX J171 3.7-3946 above 10 TeV becomes 
rather steep with photon index pa 3 (lAharonian et al.l 120071 ). 



In this Letter we suggest a search for multi-TeV 7-rays generated by the CRs that escape 
the SNR. A molec ular cloud located close to the SNR can provide an effective target for pro- 
ducti on of 7-rays ( Aharonianl Il99ll ; lAharonian &: Atoyanl Il996l ; iGabici. Aharonian &: Blasi 
20071 ). The highest energy particles (~ few PeV) escape the shell first. Moreover, generally 
they diffuse in the interstellar medium faster than low energy particles. Therefore they arrive 
first to the cloud, producing there 7-rays and neutrinos with very hard energy spectra. Note 
that an association of SN Rs with clouds is naturally expected, especially in star forming 
regions (jMontmerld 1 19791 ) . The duration of 7-ray emission in this case is determined by the 
time of propagation of CRs from the SNR to the cloud. Therefore 7-ray emission of the 
cloud lasts much longer than the emission of the SNR itself. This makes the detection of 
delayed 7-ray and neutrino signals more probable. The detection of these multi-TeV 7-rays 
from nearby clouds would thus indicate that the nearby SNR in the past was acting as an 
effective PeVatron. 



2. The model 

Consider a supernova of total energy 10 5 E51 erg exploding in a medium of density n. 
The initial shock velocity is 10 9 Wg cm/s and remains roughly constant until the mass of the 
swept up material equals the mass of the ejecta. This happens at a time ~ 200 [£51 /(n Mg)] 1 / 3 yr, 
when the shock radius is ~ 2.1 [£51/ (n Ug)] 1 / 3 pc. Then the SNR enters the Sedov phase and 
the shock radius and velocity scales with time as R s h oc t 2 / 5 and u s h oc t~ 3//5 . 
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The spectrum of particles accelerated at the SNR shock is determined by the transport 
equation: 

where D = D(p) is the momentum dependent diffusion coefficient and u the flow velocity. 
For a strong shock with compression factor r s = 4, the test particle theory predicts an 
universal shape for the CR spectrum at the shock fo(p) oc p~ A (IMalkov fc Druryll200ll ). If 
the shock is an efficient accelerator (as SNR shocks are believed to be), the CR pressure 
modifies the flow structure, making the shock more comp ressible and the spectr um of the 
accelerated particles harder, fo(p) oc p~ a with 3.5 < a < 4 (IMalkov fc Druryll200ll ). Detailed 
calculations compared with mul t iwavelength observation s of SNR suggest the values r s « 7 
and a ~ 3.7 (lEllison at al.l 120071 ; iBerezhko fc Volkil2006l ). which we adopt in the following. 



The maximum momentum of the accelerated particles is determined by a simple con- 
finement condition, namely that the diffusion length of the particles cannot exceed the 
characteristic size of the system R s h\ 



D{p r , 



<R 



U s h 



sh ■ 



(2) 



The maximum possible energies are achieved when the acceleration proceeds in the Bohm 
diffusion limit, D oc p/B s h, with B s h the magnetic field strength at the shock. In this case 
the maximum momentum decreases with time as p ma x{t) oc B s ht~ l * b . In fact, the drop of 
Pmax is even faster, given that the magnetic field is also expected to decrease with time. 
This implies that at any time, particles with momentum above p m ax{t) quickly escape the 
remnant, generating a cutoff in the s pectrum. The spectrum of t he runaway particles can 



be approximated as a 5-function (see lPtuskin fc Zirakashvilj 120051 ): 



qesc(p, t) = ~5{p - Pmax) 



d A R 



dpr, 



+ 



Vu 



-Pr, 



f (Pmax: R) 



(3) 



M, <)l 

where the integration has to be performed where the integrand is negative. Thus, to calculate 
the flux of the runaway particles one has to know: (i) the CR particle distribution function 
at Pmax at any location in the SNR, (ii) the flow velocity both inside the shock and outside 
it, where the C R precursor forms and ( in ) how the maximum momentum varies during the 



SNR evolution. iPtuskin fc Zirakashvilj (120051 ) showed that it is straigthforward to derive 



(%) and (ii) using an approximate (but reasonably accurate) linear velocity profile inside the 
tin 

ram pressure gu 2 sh and that fo(p v 



SNR and assuming that the CR pressure at the shock P^ h R is a fraction £cr of the incoming 



oc Pg R 
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The determination of p ma x and its evolution with time requires the knowledge of the 
diffusion coefficient (see Eq. [2]), which is in turn determined by the level of magnetic turbu- 
lence generated by accelerated particles themselves. This makes the problem nonlinear and 
very difficult to be solved. The value of p max depends on a few crucial but poorly known 
aspects of the problem, including the natu re of CR-driven instability operating in the shock 



precursor and the level of wave damping (IBell fc Lucekll200ll; iPtuskin fc Zirakashvifil 12005 
Blasi. Amato fc Caprioli 2007 ; Vladimirov. Ellison &: Bykov 12006 ). Because of these uncer- 
tainties, we adopt here a phenomenological approach, namely we parametrize the maximum 
momentum as p ma x(t) oc t~ 5 . We further assume p ma x ~ 5 PeV and ~ 10 GeV at the 
early and late epochs of the Sedov phase respectively. This requires 5 ~ 2.4. Remarkably, 
if the maximum momentum is a power law function of time, the spectrum of the escap- 
i ng particles integrated over the whole Sedov phase is also a power law of the form oc p~ A 
( Ptuskin Sz Zirakashvili 2005h . which is close (sligthly harder) to what needed to fit the CR 
data below the knee (IBerezinskii et al.lll990l ). 



The spectrum of the CRs inside the SNR f in (R,p,t) can be obtained from Eq. [T] by 
dropping the diffusion term, while the spectrum of the runaway CRs at a given distance R 
from the SNR and at a given time t can be obtained by solving the diffusion equation: 



df, 



out 



dt 



(R,p,t) = D IS M(p)V 2 f ou t(R,p,t) + q esc (p,t)6(R). 



(4) 



The diffusion coefficient D ISM (p) describes the propagation of CRs in the galactic disk. The 
available CR data require a power-law energy de pendence, Dtsm(E) (x E~ s , with Djsm ~ 
10 28 cm 2 /s at E ps 10 GeV and s ~ 0.3 4- 0.7 (jBerezinskii et al.lll99ol ). The constraints 
on the diffusion coefficient are obtained from the comparison between diffusion models and 
CR data and have to be considered as average galactic values. However, the conditions 
might be rather different in regions close to CR sources, in particular due to the presence 
of strong gradients in the CR distribution, whi ch may enhan ce the generation of plasma 
waves and thus suppress the diffusion coefficient (jWentzellll974l ). Below we assume D IS m = 
3 x 10 29 (E /lPeV)°' 5 cm 2 /s. The change of s within the allowed range or the choice of a 
different normalization for Djsm does not alter qualitatively the results, the main effect 
being that the characteristic time scales of the problem change prortional to 1/Djsm- 

The functions f in and f out can be used to evaluate the 7-ray and neutrino fluxes due 
to CR interactions in the ambient gas, both from the SNR itself and surrounding dense 
environments (e.g. from nearby molecular clouds). 
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3. Results 

The top panel of Fig. [1] shows the energy spectrum of 7-ray emission from interactions o 
accelerated protons with ambient medium, calculated for typical parameters characterizing 
SNRs: E 5 i = n = u 9 = 1. The bottom panel shows the emission from a cloud of mass M c i = 
1O 4 M lo cated at a distance d r j = 100 p c away from the SNR. Spectra have been calculated 



following iKelner. Aharonian fc Bugayovl (120061 ). The distance of the SNR is assumed D = 1 



kpc and different curves refer to different times after the supernova explosion. The efficiency 
of CR acceleration at the SNR shock is regulated by the parameter £cr (the ratio between 
the CR pressure at the shock to the shock ram pressure), which is assumed to be equal to 
0.3 and constant during the SNR evolution. This assumption is rea sonable for strong shocks 



for w hich the acceleration efficiency saturates to very high values (IBlasi. Gabici fc Vannoni 



20051 ) . and it becomes less reliable in the late stages of the Sedov phase, when the SNR shock 



becomes progressively weaker. 

Early in the Sedov phase (curve 1, 400 yr after the explosion), the 7-ray spectrum from 
the SNR is hard and extends up to > 100 TeV, revealing the acceleration of PeV particles. 
The hardness of the spectrum reflects the fact that, due to nonlinear effects in particle 
acceleration, the underlying CR spectrum becomes harder than p 2 fo(p) oc p~ 2 . Conversely, 
the 7-ray flux from the cloud is extremely weak, because for the epoch of 400 yr after the 
explosion CRs do not have sufficient time to reach the cloud. The emission of > lOOTeV 
photons from the SNR lasts a few hundreds years, and after that the cutoff in the 7-ray 
spectrum moves to lower energies (curves 2, 3 and 4 correspond to the epochs of 2 10 3 , 8 10 3 
and 3.2 10 4 yr after the explosion). As time passes, CRs finally reach the cloud and produce 
there 7-rays when interacting with the dense cloud environment. This makes the cloud an 
effective multi-TeV 7-ray emitter, with a flux at the sensitivity level of next generation 
Cherenkov telescopes operating in that energy range. As lower and lower energy particles 
reach the cloud, the peak of the 7-ray emission accordingly shifts towards the lower, TeV and 
GeV, energies at flux levels which can be probed by ground based instruments and GLAST. 

The shape of the 7-ray spectrum is naturally explained as follows: at a time t, only 
particles with energy above E*, given by d c i m ^yGD IS Ai(E*)t, reach the cloud. Thus the 
CR spectrum inside the cloud has a sharp low energy cutoff at E*. The corresponding 7-ray 
spectrum peaks at the energy p» 0.1E* and the spectral slopes above and below the peak are 
~ E~ 2 ~ s and ~ E^ 1 respectively. 

The multi-TeV hadronic 7-ray emission from the cloud is significantly weaker than the 
one from the SNR, but its detection might be easier because of its longer duration (< 10 4 
yr versus few hundreds years). Moreover, the leptonic contribution to the cloud emission is 
likely to be negligible. Electrons accelerated at the SNR cannot reach the cloud because they 
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remain confined in the SNR due to severe synchrotron losses. Secondary electrons can be 
produced in the clou d, but they cool mainly via synchrotron emission in the cloud megnatic 



field ~ 10 -j- lOOyuG (jCrutcherl Il999l ). This makes the production of > TeV 7-rays due to 



inverse Compton scattering and non-thermal Bremsstrahlung negligible. 

Fig. [2] shows the muonic neutrino fluxes from the SNR (top panel) and the cloud 
(bottom panel) for the same parameters adopted in Fig. HJ The flux at Earth is a factor of 
Ri 2 smaller than what showed due to neutrino oscillations. For a young SNR (curve 1), the 
spectrum extends up to ~ 100 TeV, where km 3 -scale neutrino telescopes achieve their best 
performance. The flux level is > 10~ 11 TeVcm -2 s~ 1 , which makes such sources detectable 
in a few years. Unfortunately, the high energy cutoff in the neutrino spectrum moves fast 
towards low energies (curves 2 and 3), making the detection more problematic. Note that the 
7-ray and neutrino fluxes from the SNR itself and from the cloud scale as n/D 2 and M c i/D 2 
respectively. Thus, the presence of a very massive molecular cloud with mass M c \ ~ 1O 5 M 
would considerably increase the chances of detection of multi-TeV 7-rays and neutrinos from 
such systems. 

The impact of the distance between the SNR and the cloud on the flux of 7-ray and 
neutrino emission is illustrated in Fig. [3X where 7-ray (thick lines) and neutrino (thin lines) 
spectra are shown for a cloud at a distance of 50 and 200 pc from the SNR. Note that the 
neutrino flux of the cloud may become marginally detectable only if the cloud is very close 
to the SNR. On the other hand, clouds can be detectable in 7-rays even if their distance 
from the SNR is as large as ~ 200 pc (right panel). Remarkably, the TeV emission from the 
cloud lasts > 10 4 yr, significantly longer than the emission from the SNR. 



4. Conclusions 

The acceleration of CRs up to the knee in SNRs can be revealed by means of observations 
of mufti- TeV 7-rays and neutrinos coming from the SNR and nearby molecular clouds. The 
emission from the clouds is weaker than the one from the SNR, but lasts longer, enhancing 
the probability of detection. Both 7-rays and neutrinos are emitted with fluxes detectable 
by currently operating and forthcoming instruments. Since the 7-ray spectra from clouds 
are extremely hard, 7-ray telescopes operating at very high energies (> 10 TeV) would be 
the best instruments for this kind of study. A detection of such emission would both reveal 
the acceleration of PeV CRs and suggest the best targets for neutrino observations. 



Finally, we note that Yamazaki et all 1 
unidentified TeV sources (lAharonian et al 



(20061) recently proposed an asso ciation between 
J 120021 l2006bl : lAbdo et all l2007i ) and old SNRs. 
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Such idea seems to us questionable due to the fact that old SNRs cannot confine multi-TeV 
particles and thus emit effectively 7-rays only at sub-TeV energies (see Fig. [1]). In the 
scenario we proposed in this paper, unidentified TeV sources might still be associated with 
old SNRs, the gamma ray emission being produced in the interaction of escaping CRs with 
nearby dense clouds. 
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Fig. 1. — Gamma ray spectra from the SNR (TOP) and from a cloud of 1O 4 M located 100 
pc away from the SNR (BOTTOM). The distance is 1 kpc. Curves refer to different times 
after the explosion: 400 (curve 1), 2000 (2), 8000 (3), 3.2 10 4 (4) yr. 
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Fig. 2. — Muonic neutrino spectra from the SNR (TOP) and from the cloud (BOTTOM) 
for the same set of parameters of Fig. [TJ Neutrino oscillations are not taken into account. 
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Fig. 3. — Gamma ray (thick lines) and neutrino (thin lines) spectra from a cloud located at 
50 (left) and 200 pc (right) from the SNR. Different curves refer to different times after the 
explosion. 



